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ABSTRACT 

We present a statistical comparison of three different estimates of cluster 
mass, namely, the dynamical masses obtained from the velocity dispersion of 
optical galaxies, the X-ray masses measured from the temperature of X-ray 
emitting gas under the assumption of isothermal hydrostatic equilibrium, and 
the gravitational lcnsing masses derived from the strong/weak distortions of 
background galaxy images. Using a sample of 29 lensing clusters available in 
literature, we have shown that the dynamical masses are in agreement with 
the gravitational lensing masses, while the X-ray method has systematically 
underestimated cluster masses by a factor 2-3 as compared with the others. 
These results imply that galaxies indeed trace the gravitational potential of 
their clusters, and there is no bias between the velocities of the dark matter 
particles and the galaxies in clusters. The X-ray cluster mass discrepancy is 
probably from the simplification in the models for the X-ray gas distribution 
and dynamical evolution. 

Subject headings: cosmology: theory — galaxies: clusters: general - 
gravitational lensing 



1. Introduction 


Clusters of galaxies are the largest coherent and gravitationally bounded objects in 
the universe. The precise determination of their gravitational masses is crucial for our 
understanding of formation and evolution of cosmic structures, for mapping of matter 
distribution on large-scales and also for measurement of the present mean mass density of 
the universe (hlo)- Historically, cluster masses are derived from the dynamical analysis of the 
observed velocity dispersion of cluster galaxies based on the virial theorem, which results 
in the so-called virial cluster masses M v . With the development of the X-ray astronomical 
techniques clusters can be selected from the X-ray emission of the hot diffuse intracluster 
gas, giving rise to the X-ray cluster masses M x when combined with the assumption 
of hydrostatic equilibrium. Over the past decade the detection of the gravitationally 
distorted images of faint distant galaxies behind some clusters of galaxies provides another 
independent mass estimate: the gravitational cluster masses Mi ens . In particular, Mi ens 
are obtained regardless of the cluster matter state and components. These three methods 
should be incorporate, and comparisons of their results would yield a very useful clue 
regarding the dynamical evolution of clusters and a test for the accuracy of cluster mass 
determinations. 

The early studies based on a few selected clusters in which both lensing and optical 
and/or X-ray data are available claimed a cluster mass discrepancy by a factor of typically 
2 ~ 3 among the three methods (Wu, 1994; Fahlman et al. 1994; Miralda-Escude & Loeb 
1995), while a consistency between dynamical masses and gravitational lensing derived ones 
has been also reported in some cases, e.g. PKS0745-191 (Allen, Fabian & Kneib 1996). 
We have recently carried out a statistical comparison of the overall cluster radial matter 
distributions determined from X-ray observation and gravitational lensing and concluded 
that the X-ray analysis may have systematically underestimated cluster masses at least in 
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the central regions (Wu & Fang 1996; hereafter Paper I). 

From the theoretical point of view, the above three methods have to give the same 
cluster masses if clusters are dynamically relaxed. This probably accounts for the result 
in clusters like PKS0745-191 which appear to be regular in optical/X-ray morphologies. 
Though one may attribute the reported mass discrepancy to the nonthcrmal pressure (Loeb 
& Mao 1994; Ensslin et al. 1997) or the projection effect (Miralda-Escude & Babul 1995; 
Cen 1996), it is most likely that the problem is relevant to the cluster matter distribution 
and dynamical evolution. Indeed, both optical and X-ray observations have revealed the 
presence of substructures in most of clusters, implying that clusters may be still in the era 
of formation. In particular, the recent spatially-resolved measurements of gas temperature 
in some clusters illustrate the complex two-dimensional patterns including the asymmetric 
variations and the significant decline with radius (e.g. Henriksen & White 1996; Henriksen 
& Markevitch 1996; Markevitch 1996), which are the strong indicators of the effect of 
substructure merging. For those clusters, we are unable to apply the equation of hydrostatic 
equilibrium to the X-ray emitting gas. Therefore, it is naturally expected that the X-ray 
cluster mass obtained under the assumption of a hydrostatic equilibrium is in principle 
not representative of the true cluster mass, and thereby should be different from the 
gravitational lensing derived mass and/or even the virial mass. This scenario is supported 
by the recent numerical study of X-ray and lensing properties of clusters of galaxies from 
the standard CDM simulations by Bartelmann & Steinmetz (1996). They found that the 
cluster masses can be biased low if the traditional /3 model is adopted for the intracluster 
gas, which is especially true when clusters exhibit pronounced substructures, while the 
strong lensing preferentially selects the clusters that are dynamically more active than the 
average. 

Alternatively, it should be noticed that M v derived from the virial theorem and M x 
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derived from the hydrostatic equilibrium have very different physical implications. M v is 
related to the galaxy velocity dispersion and galaxy number density. Comparison of M v 
and Mi ens would set constraints on the bias parameter between the velocities of the galaxies 
and of the dark matter and test whether optical galaxies trace the gravitational potential of 
the cluster. On the other hand, M x depends on the temperature variation and the density 
profile of the hot diffuse gas, which may suffer from the influence of the possible existence 
of the turbulence and magnetic field. As a result, comparison of M x and M v may allow 
one to determine how significant the nonthcrmal pressure contributes to the computation 
of M x because galaxies are unaffected by the nonthermal pressure in clusters. Furthermore, 
the comparison of the three cluster estimates may help to solve the puzzle of the baryonic 
matter excess in clusters if clusters provide a fair sample of the universal baryon fraction 
(see Paper I for summary). Recall that the baryon fractions in clusters are computed using 
the dynamical masses M x and/or M v , which may have large uncertainties if most of clusters 
are still in the process of violent merging. As a whole, galaxies and gas particles have 
probably experienced very different evolutions in the formation of clusters and then exhibit 
different dynamical states and density distributions in clusters today. 

Numerous lensing and optical/X-ray observations of clusters have now made it timely 
and possible to carry out these comparisons. Unlike Paper I that chose X-ray and lensing 
data separately from literature, we now work with the lensing clusters only, in which the 
strongly and/or weakly distorted images of background galaxies have been observed. While 
this paper was in the refereeing stage, we received a preprint by Small et al. (1997) who 
made a similar but sophisticated investigation for 11 distant clusters observed with HST. 
Throughout this paper we adopt a matter-dominated flat cosmological model of Oo = 1 and 
a Hubble constant of Hq = 50 km s^ 1 Mpc. 
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2. Cluster sample and mass estimates 

Strongly/weakly distorted images of background galaxies have been so far detected 
in ~ 40 clusters of galaxies (Fort & Mellier 1994; Wu 1996; Van Waerbeke & 

Mellier 1996). For our purpose we only select those lensing clusters whose X-ray 
luminosity(L, x )/teniperature(T) and/or optical galaxy velocity dispersion (a) are available 
in literature. This results in a non-exhaustive list of lensing cluster sample (Table 1) 
containing 29 clusters and 39 measurements. 

EDITOR: PLACE TABLE | HERE. 

Temperature of the hot diffuse gas is the most important parameter for the X-ray 
cluster mass determination. Unfortunately, only few clusters in Table 1 have the measured 
temperature: A370 (T = 8.8 ± 0.8 keV), A1689 (T = lO.l/j/g keV), A2163 (T = 13.9 keV), 
A2218 (T = 6.7 keV), MS0451(T = 10.4 ± 1.2 keV) and PKS0745 (T = 8.6keV). This 
arises from the difficulty of the X-ray spectroscopic observations. Nonetheless, a correlation 
between cluster temperature and bolometric X-ray luminosity L Xjbo io has been well 
established (e.g. Edge & Stewart 1991; Fabian et al. 1994): L X)bo i 0 = io 43 - 06±0 - 08 T 2 ' 68±010 
while L x in unit of erg s _1 can be relatively easily obtained. Such a relationship thus 
enables us to translate L x into T. Assuming a mechanism of free-free bremsstrahlung for 
the gas-emission and adopting an approximate Gaunt factor given by Mewe et al. (1986), 
we have computed the cluster temperature T for each cluster in Table 1. Note that the 
error bars in the resulting T account for both the uncertainties of the L x measurement and 
of the L x -T relation. It is seen that the agreement between the measured temperature and 
the estimated one is fairly good for the above six clusters. 

Cluster optical and X-ray morphologies reflect their dynamical state at different 
evolutionary stages, which turn to be quite varied observationally. Apparently, a well- 
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relaxed system would appear to be regular. MS0440 and PKS0745 in Table 1 might be 
two good examples of this kind. On the other hand, clusters that are currently forming 
may have very irregular morphologies. Substructures in optical/X-ray images are the most 
common feature in those clusters and on-going collisions of the substructures are the most 
convincing evidence for interaction. Another indicator for the cluster non-hydrostatic 
process is the complex temperature patterns observed recently from a number of spatially 
resolved measurements of the gas temperature. All these sorts of clusters are marked by “I” 
in Table 1. Moreover, we utilize “E” to denote the rest clusters in which substructures do 
not clearly present but their X-ray/optical maps look more or less like ellipses in shape. 


Fitting the azimuthally averaged X-ray surface brightness profile by the usual (3 model 
and assuming a hydrostatic equilibrium for the X-ray emitting gas, we can obtain the 
projected X-ray cluster mass m x within cluster radius r through (Wu 1994) 


m x = 1.13 x 10 u Mo/3f it 


kT 


Mpc / \ keV, 


m[r 


( 1 ) 
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and r xc and R are the core radius in /3 model and the cluster physical radius, respectively. 
A straightforward computation shows that m depends very weakly on R and we will thus 
take a value of R = 3 Mpc in the following calculation. 


We now give another way to estimate the cluster mass. We model the cluster matter 
distribution by an isothermal sphere which is characterize by a core radius r<j c and the 
optical galaxy velocity dispersion er. We use the term “virial” or “dynamical” cluster mass 
to denote the mass given by the optical galaxy velocity dispersion, though this differs from 
the usual virial analysis which utilizes the galaxy number density in cluster rather than 
the dark matter profile. Our attempt here is to examine whether a can provide a good 
mass estimate when compared with gravitational lensing method. Yet, the radial galaxy 



distribution in cluster is not well constrained. Recall that the traditional King model has 
been questioned (Bahcall & Lubin 1994). For a softened isothermal sphere, the projected 
“virial” mass within a radius of r is simply 


rn v 


7T <7 


G 



( 2 ) 


Finally, with gravitational lensing method one is able to determine the projected 
gravitational cluster mass within the arc position r arc or the distance r from the cluster 
center on the cluster plane: 

{ n r lrc^crit, for arc/arclet; 

(3) 

7rr 2 (j(r)E c „ t , for weak lensing, 

where Y> crit = (c 2 / AnG)(D s /D^Dds) is the critical surface mass density with Da, D s and 
Dds being the angular diameter distances to the cluster, to the background galaxy and from 
the cluster to the galaxy, respectively, and £( r ) measures the statistics of the shear held 
7 of background galaxy images induced by the cluster (Fahlman et al. 1994). Note that 
in eq.(3) arclike images have been presumed to trace the Einstein radius of the cluster so 
that the alignment parameter is approximately taken to be zero. Furthermore, we need to 
make two remarks: (l)Redshift data are still not available for nearly half of the arcs listed 
in Table 1, for which we have assumed a redshift of z s = 0.8. This leads to the decrease of 
mass estimate mi ens by a factor of 1.4 for a typical arc-cluster at 27 = 0.3 if the background 
galaxy is set to be z s = 2. (2)Weaking lensing analysis provides only a low bound on mi ens . 

The resulting three mass estimates have been given in Table 1. Uncertainties in most of 
the mi ens results and some of the m v results are hard to evaluate at present. Alternatively, 
in the computation of rn x we adopt a mean core radius of r xc = 0.25 Mpc but allow r xc to 
vary from 0.1 Mpc to 0.5 Mpc. So, the error bars in the final result of m x include both the 
uncertainties in T and the ones in r xc . Moreover, all the results of m x in Table 1 correspond 



to (3fit = 1 while the observationally fitted value is around /3f it ~ 0.67 (e.g. Jones & Forman 
1984). For the dynamical cluster mass rn v given by the optical galaxy velocity dispersion, 
we have assumed rdc = 0. 


3. Results and discussion 

Comparisons of cluster mass estimates from X-ray gas hydrostatic equilibrium, 
dynamical analysis and gravitational lensing are shown in Figure 1(a) and (b) using the 
data of Table 1. An immediate conclusion is that mi ens agrees essentially with m v while a 
systematic excess of mi ens with respect to rn x is detected. This can be clearly demonstrated 
by the following best-fit relation to the data: 

miens = (1-42 ± 0.99)771, = (2.23 ± 1.15 )/3j£m x , (4) 

in which the uncertainties are the lcr errorbars. Meanwhile, Fig. 1(a) and (b) also illustrate 
the influence of cluster morphologies on the relations between mi ens and m v and between 
mi ens an d m x , respectively. Apparently, it is very unlikely that cluster morphologies can 
lead to a remarkable difference in the results. 

EDITOR: PLACE FIGURE 0 HERE. 

In Fig.2(a) and (b) we display the variations of mi ens , m v and m x with the cluster 
radius. Basically, Fig.2(b) provides a result similar to the one of Paper I, i.e., there 
is a systematic discrepancy between m x and mi ens inside the cluster core radius, and 
the projected gravitational cluster mass obtained with lensing mothed follows a simple 
power-law of ~ r 1 ' 3 . Again, in contrast to m x , m v agrees statistically with mi ens over all the 
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EDITOR: PLACE FIGURE @ HERE. 

Additionally, we have computed the [3 parameter characterizing the specific energies of 
the galaxies and the gas in clusters, (3 spec = cr 2 / (kT /fjLm p ) where /im p = 0.59 is the mean 
particle mass. Our best-fit value with the 18 data points (both T and cr are available) in 
Table 1 reads (3 spec = 1.29 ± 0.71, while f3 spec reduces to 1.17 ± 0.50 if AC114 is excluded. 
The best-fit relation between the galaxy velocity dispersion and the gas temperature is 

(cr/km s" 1 ) = 10 2 - 64±0 ' 11 (T/keV)°- 51±0 - 13 . (5) 

It appears that our best-fit average f3 spec and cr-T relation are consistent with the previous 
work [see Girardi et al. (1996) for summary]. Based on such a good fitness of eq.(5) alone, 
one might conclude that the galaxies and the gas are in hydrostatic equilibrium with the 
same cluster potential, as was claimed by Lubin & Bahcall (1993). However, our result of 
eq.(4) raises a new question of how one could reconcile the discrepancy between mi ens (or 
m v ) and m x with the good correlation between a and T. 

Alternatively, it turns out from Fig. 1(a) that a zero core radius r dc = 0 for the 
cluster mass profile provides a good fit to the lensing data. The best-fit core radius r dc by 
requiring m v = mi ens is r dc = —0.09 ±0.24 Mpc, indicative of rather a compact dark matter 
distribution in clusters. This result is compatible with the early studies of giant arcs and 
statistical lensing of arcs/arclets which report a small core radius (< 0.1 Mpc) for the dark 
matter profile of the arc clusters (e.g. Hammer 1991; Wu & Hammer 1993; Grossman & 
Saha 1994). 

While there is a significant evolution of X-ray luminosity clusters with redshift (Edge 
et al. 1990; Gioia et al. 1990; Henry et al. 1992), the deficit of the X-ray cluster mass may 
be relevant to the cluster evolution. Recall that the local L x -T relation established at low 
redshift z < 0.1 was employed to estimate the cluster temperature, whereas most of the 
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clusters in Table 1 are actually located at intermediate redshift z ~ 0.2-0.5. Therefore, 
it would be useful to examine the dependence of the ratio of rrii ens to rn x on the cluster 
redshift. Our best-fit relation is 

m lens _ ^g0.09±0.10/-^ _|_ „,)1.7±0.9 (g) 

m x 

Namely, the cluster mass discrepancy between the X-ray analysis and the gravitational 
lensing method is indeed related to the evolutionary history of clusters. Since m x is 
proportional to T according to eq.(l), the cluster temperature has a similar variation with 
redshift. This yields a temperature ratio of 1.5i(j;§ f° r cluster at redshift z = 0.33 relative to 
the one at z = 0.035, in consistent with the result (1.41 q; 3 ) of Henry et al. (1994). Though 
the cluster temperature evolution since intermediate redshift is moderate, it may account 
for the mass discrepancy we report in the present paper. 

We now discuss briefly the significance of the consistency/discrepancy between 
the cluster masses derived from gravitational lensing, dynamical analysis and X-ray 
observations. 

Both gravitational lensing and “virial” methods yield nearly the same cluster masses, 
which have several implications: First, galaxy velocity dispersion indeed provides a 
good estimate of cluster mass. Second, the dark matter particles and the galaxies have 
approximately the same velocity dispersion, i.e., there is no velocity bias in clusters of 
galaxies. Third, mass follows the light. These arguments are comparable with the recent 
dynamical analysis of the CNOC cluster sample (Carlberg et al. 1996), which has found 
strong evidence that galaxies are effectively in equilibrium with their host cluster. However, 
our finding disagrees with the numerical result that the velocity biasing parameter is 
~ 0.7-0.8 (Carlberg & Dubinski 1991; Couchman & Carlberg 1992). 

On the other hands, cluster mass estimate based on the X-ray temperature assuming 
an isothermal hydrostatic equilibrium has systematically underestimated cluster masses 
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by a factor of ~ 2-3 as compared with gravitational lcnsing and “virial” method, which 
demonstrates that the gas particles may not be a good tracer of the gravitational potential 
of the cluster. The recent high-spectral resolution observations do reveal the occurrence of 
complex temperature maps, indicating the on-going merger activities in clusters and the 
cluster evolution with cosmic epoch. Overall, the simplification of modeling the temperature 
(isothermal) and gas density profile (spherical) is responsible for the deficit of the X-ray 
cluster mass detected in this paper. One may also attribute the mass disprepancy to the 
nonthermal pressure in clusters (Loeb & Mao 1994; Ensslin et al. 1997), which affects 
the gas particles while produces no effect on galaxy distribution and velocity dispersion. 
Finally, it is pointed out that our results still contain large scatters due to the scarcity of 
the lensing data, and a large cluster sample will be needed to confirm our finding. 

Useful discussion with Frangois Hammer and Shude Mao and valuable suggestions by 
an anonymous referee are gratefully acknowledged. WXP wishes to thank the hospitality 
of DAEC, Observatoire de Paris-Meudon, where part of this research was made. This work 
was supported by the National Science Foundation of China. 
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Fig. 1.— Gravitational lensing derived cluster mass mi ens plotted against dynamical cluster 
mass rn v (a) and X-ray cluster mass m x (b). Additionally, cluster morphologies are 
represented by different symbols: open circles - regular, filled circles - elliptical and asterisks 
- irregular. The dashed lines are the least square fits of the data to a power-low. 

Fig. 2.— Comparisons of cluster radial mass distributions derived from gravitational 
lensing and dynamical method (a) and from gravitational lensing and X-ray gas hydrostatic 
equilibrium (b). Note that the (3 paramter is taken to be [3f lt = 0.67 in (b). 
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Table 1: Lensing cluster sample 
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I 

1 6 +0 ' 5 

4 ' D —0.8 

0.825 

arc 

0.93 

0.14 

1.6 

12 

(jiiiqooo 

0.316 

(l-)(0.1 — 2.4) 

(2.8±g;*) 

1375 

I 

(0.38+ s lg ) 

2.0 

arc 

0.913 

0.15 

1.9 

13 

(00024 

0.391 

2-7(2-io) 

4 7+0- 3 
^•'-0.2 

1300 

E 

-I 9+0.5 

1-Z —0.5 

2.8 

arc 

1.39(?) 

0.22 

3.6 

14 

-fa 

n 

<gg0302 







37 

w.l. 


3.0 

40 

15 

0.423 

4-8(0.7—5) 

6 O” 1 " 0,3 
°‘ u -0.3 


I 

0 5 fi+°- 49 
U.OD_o 28 


arc 


0.12 

1.6 

16 

pt2244 

0.328 

1-5(0.5-4.5) 

4 4+ 0 - 2 


I 

0 11+ 015 
U-1± -0.06 


arc 

2.237 

0.061 

0.25 

17 

^«0440 

0.190 

4-01(0.3-3.5) 

4 9+°- 8 
*- y -0.6 


R 

0 26 +0 ‘ 34 
u - zd - 0.14 


arc 

0.53 

0.089 

0.89 

18 

g0451 

0.55 

20(0.3-3.5) 

9 7+ 2,0 
y ' ‘ —1.5 

1371 

I 

1 9+ 1 ' 4 
1 - y -i.o 

2.61 

arc 


0.19 

5.2 

19 

M31006 

0.221 

4-82(0.3_3.5) 

E q + 0.9 
— 0.7 

906 

I 

n 99+O.32 
U-ZZ —0.12 

0.47 

arc 


0.079 

0.57 

19 







0 63 +0 ‘ 60 

U.OO_o 33 

0.84 

arc 


0.14 

1.8 

19 







1 9+0- 9 
±,y —0.9 

1.68 

arc 


0.28 

7.2 

19 

MS1008 

0.301 

4-49(0.3-3.5) 

c 1+0.9 
°- 1 -0.7 

1054 

I 

1 7+O.9 
' —0.8 

2.11 

arc 


0.26 

6.1 

19 

MS1054 

0.826 

9-02(0.3-3.5) 

6 q + 1 -3 


I 

20.211° 


w.l. 


1.9 

28+ 

20 

MS1224 

0.327 

4-61(0.3-3.5) 

E rj + 0.9 
°' z -0.7 

802 

E/I 

8-1+7 

4.5 

w.l. 


0.96 

7.0 

21 

MS1445 

0.259 

16.0(0.3-3.5) 

8.8+;® 

1133 

I 

0-55±“;” 

0.92 

arc 


0.098 

0.86 

19 

MS1621 

0.426 

4-55(0.3-3.5) 

5-i+o: 9 7 

793 

R 

0 076”*~ 0 ' 13 
u.u/D _ 0 043 

0.21 

arc 


0.046 

0.24 

22 

MS1910 

0.246 

4.38(0.3-3.5) 

c n+0.9 
O,U -0.7 


I 

Q q + 0.9 
z, °— 0.9 


arc 


0.33 

9.6 

19 

MS2053 

0.583 

5-8(0.3-3.5) 

s+o.s 


R 

0 51 + 0- 57 
U - O 1 - 0.28 


arc 


0.119 

2.6 

22 

MS2137 

0.313 

15-6(0.3-3.5) 

8 - 7+3 

960 

R/E 

0 44+ 0 - 52 
u, **-0.25 

0.59 

arc 


0.088 

0.71 

23 

MS2318 

0.187 

6-84(0.3_3.5) 

6 1+ 11 

1 — 0.9 


I 

0 515+O 6I 
u.OQ-o 30 


arc 


0.12 

1.3 

19 

AC114 

0.31 

4-0(0.1-2.4) 

s-otg:? 

1649^220 

I(?) 

2 5 +1 -° 
Z *°-1.0 

6.9+° 

arc 

0.639 

0.35 

13 

24 

PKS0745 

0.103 

28(2-10) 

10.6tg;? 


R 

n Q 9 + 0 .O 8 

u -°^— 0.05 


arc 

0.433 

0.046 

0.30 

25 

RXJ1347 

0.451 

62(0.1-2.4) 

17-5^'j 


I 

c q+ 3.8 
°- u -2.6 


arc 


0.24 

6.6 

26 






-14- 


°X-ray luminosity in unit of 10 44 erg s _1 measured in the e energy band (keV). 
Estimated temperature (keV) from the L x -T relation. 

°Cluster X-ray/optical morphology: E-elliptical, I-irregular and R regular. 
^Projected cluster mass within radius r in unit of 10 14 Mq. 
e Lensing phenomenon: arc(s) or weak lensing (w.l.). 

■^Distance from the center of cluster where arc(s) and/or weak lensing are detected. 
9 Only one reference is provided for lensing features. 
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al. 1994; (22)Luppino & Gioia 1992; (23)Mellier et al. 1993; (24)Smail et al. 1995a; (25)Allen et al. 1996; 
(26)Schindler et al. 1995. 



